Quantum-limited shot noise in graphene 
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We calculate the mode-dependent transmission probability of massless Dirac fermions through an 
ideal strip of graphene (length L, width W, no impurities or defects), to obtain the conductance 
and shot noise as a function of Fermi energy. We find that the minimum conductivity of order 

/h at the Dirac point (when the electron and hole excitations are degenerate) is associated with 
a maximum of the Fano factor (the ratio of noise power and mean current). For short and wide 
graphene strips the Fano factor at the Dirac point equals 1/3, three times smaller than for a Poisson 
process. This is the same value as for a disordered metal, which is remarkable since the classical 
dynamics of the Dirac fermions is ballistic. 

PACS numbers: 73.50.Td, 73.23.-b, 73.23.Ad, 73.63.-b 



Two recent experiments 0, 01 have discovered that the 
conductivity of graphene (a single atomic layer of carbon) 
tends to a minimum value of the order of the quantum 
unit e^//i when the concentration of charge carriers tends 
to zero. This quantum-hmited conductivity is an intrinsic 
property of two-dimensional Dirac fermions (massless ex- 
citations governed by a relativistic wave equation), which 
persists in an ideal crystal without any impurities or lat- 
tice defects H, IM ly| ■ In the absence of impurity scat- 
tering, and at zero temperature, one might expect the 
electrical current to be noiseless. In contrast, we show 
that the minimum in the conductivity is associated with 
a maximum in the Fano factor (the ratio of noise power 
and mean current). The Fano factor at zero carrier con- 
centration takes on the universal value 1/3 for a short 
and wide graphene strip. This is three times smaller than 
the Poissonian noise in a tunnel junction and identical to 
the value in a disordered metal 7, 8] — even though the 
classical dynamics in the graphene strip is ballistic. 

Shot noise measurements have proven to be a valu- 
able diagnostic tool in carbon nanotubes, which can be 
thought of as roUed-up sheets of graphene. Very low 
shot noise in well-contacted bundles of single-wall nano- 
tubes is an indication of nearly ballistic one-dimensional 
transport 9]. Super-Poissonian noise has been found 
in a quantum dot formed out of a single-wall nano- 
tube, and explained in terms of inelastic tunneling in 
this zero-dimensional system [Tol | . Our prediction of sub- 
Poissonian shot noise in two-dimensional graphene is an- 
other manifestation of the importance of dimensionality 
for quantum transport. 

Our analysis of the shot noise was inspired by an in- 
sightful recent paper of Katsnelson who used the 
Landauer transmission formula to obtain the quantum- 
limited conductivity. Following the same approach, we 




*On leave from: Marian Smoluchowski Institute of Physics, Jagiel- 
lonian University, Krakow, Poland. 



FIG. 1: Schematic of a strip of graphene of width W , con- 
tacted by two electrodes (black rectangles) at a distance L. A 
voltage source drives a current through the strip. A separate 
gate electrode (not shown) allows the carrier concentration in 
the strip to be tuned around the neutrality point. 



calculate the transmission probabilities of Dirac parti- 
cles through a strip of graphene in the geometry of Fig. 
n (An earlier study of the same geometry counted the 
number of propagating modes, without determining their 
transmission probabilities |llj|.) The result depends on 
the aspect ratio W/L of the strip and also on microscopic 
details of the upper and lower edge. For short and wide 
strips [W/L ^ 1) these microscopic details become in- 
significant. For that reason we first discuss the simplest 
case of an edge which is smooth on the scale of the lattice 
spacing. This corresponds to confinement of the carriers 
by lattice straining. The opposite case of an abrupt edge 
(corresponding to confinement by etching) is considered 
later on. 

The band structure of graphene has two valleys, which 
are decoupled in the case of a smooth edge. In a given val- 
ley the excitations have a two-component envelope wave 
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function ^ = (^'i, ^'2), varying on scales large compared 
to the lattice spacing. This continuum description of the 
electronic states in graphene has been found to be quite 
accurate '12^, and we will test it later on by comparing 
with a numerical solution of the scattering in a tight- 
binding model. The two components of refer to the 
two sublattices in the two-dimensional honeycomb lat- 
tice of carbon atoms. (The additional spin degeneracy 
of the excitations does not play a role here.) The wave 
equation for ^' is the Dirac equation, 

[vp^a^j:+vpyay+v^M{y)<7,+eV{x)]'i'{r) = e*(r), (1) 

with V the velocity of the massless excitations of charge 
e and energy e, p — —ihd/dr the momentum operator, 
r — (x, y) the position, and ui a Pauli matrix. We choose 
the zero of energy such that the Fermi level is at e = 0. 

The mass term M{y) is zero in the interior of the strip 
and rises to 00 at the edges y = and y — W, thereby 
confini ng t he particles. As shown by Berry and Mon- 
dragon io\ , infinite mass confinement corresponds to the 
boundary condition 
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As a result of this boundary condition, the transversal 
momenta are quantized as 
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n = 0,1,2, 



(3) 



with n labeling the modes. The quantization condition 
for Dirac particles confined by an infinite mass differs 
from the one for normal electrons confined by an infi- 
nite potential by the offset of 1/2, originating from the tt 
phase shift in the boundary condition 

The electrostatic potential V{x) = Igatc for < a; < L, 
varied by a gate voltage, determines the concentration of 
the carriers in the strip. The value V^ato = corresponds 
to charge neutrality, being the point where electron and 
hole excitations are degenerate (known as the Dirac 
point). We model the electrodes by taking a large value 
V[x) = Voo in the leads a; < and x > L. (The parame- 
ter Voo will drop out of the results, if \Voo\ 3> |V^atc|-) 

We calculate the transmission probabilities at the 
Fermi level by matching modes at x = and x = L. 
The matching condition for the Dirac equation is the 
continuity of the two components of This ensures 
the conservation of the local current density j{r) = 
ew^'t(r) • cr • ^'(r), with cr — {ax,(Jy), without requiring 
continuity of derivatives (as needed for the Schrodinger 
equation). There is a separate transmission probability 
Tn for each of the N propagating modes in the leads, 
because the matching condition does not mix the modes. 
(The integer iV > 1 is given by = Int {kooW/n + i), 
with e|yoo| = fivkoc-) 

Details of the calculation are given in the Appendix. 
At the Dirac point V^atc = the transmission probability 



reads 

Tn = 



1 



cosh L(7„ + (gf„/fcoo)^ sinh Lg„ 

for A^ > W/L. (4) 



cosh^[7r(7i + 1/2)L/Ty] 

The formula (HJ is essentially different from the textbook 
formula 0| for the transmission probability of nonrela- 
tivistic electrons through a potential barrier, which van- 
ishes in the limit A^ — > 00 at zero energy (relative to the 
top of the barrier). 

The finite transmission probability at the Dirac point 
tends to the ballistic limit T„ — > 1 with increasing | Viatel ■ 
For N 00 we find the expression 
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3^"^{qn - kn + + 6 '^"^ {qn - kn - i^Y 



(5) 

with K — e|Vgatc|/^J^^ and kn = \J q^ — for (/„ > k or 
kn = i\/ - for q„ < k. 

The conductance G and Fano factor F follow by sum- 
ming over the modes, 



G = .go 



ji=0 



F ■ 



X^rt^O ^ra(l " Tn) 
Z^„=0 ^" 



(6) 



with go = Ae'^/h. (The factor 4 accounts for the spin and 
valley degeneracy.) The dependence of the conductivity 
a = G X L/W and the Fano factor at V^atc = on the 
aspect ratio W/L is plotted in Fig.|21(sohd curves). The 
dependence on Igatc at a fixed value of W/L is shown in 
Fig. El 

Fig. 121 also contains results for a boundary condition 
corresponding to an abrupt edge (dashed curves). We 
considered a "metallic armchair" edge, in which the car- 
bon lattice contains a multiple of three hexagons in the 
transverse direction, terminated at y = and y = W 
by a horizontal bond. This edge mixes the valleys, so 
we need to consider a four-component wave function 
^ = (*!, *2)- The first two components satisfy 

the Dirac equation (Q, without the mass term, and the 
second two components satisfy the same equation with 
Py — > —py The boundary condition is [T^ 



*i ^ ^2 = *2> at 2/ = 0, W. 



(7) 



The valley degeneracy is broken for the lowest mode 
(n = 0), which is nondegenerate, while all higher modes 
(n = 1, 2, . . .) retain the twofold valley degeneracy (over 
and above the twofold spin degeneracy, common to all 



modes). For Vg 



gate 



and A^ — > oo the transmission 



probabilities are given by 



Tn 
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cosh'' {irnL/W) 



n = 0,1,2, 



(8) 



The essential difference with the result Q for the smooth 
edge is due to the absence of the 1/2 offset in the quanti- 
zation condition of the transverse momentum. The differ- 
ent boundary condition changes the strip from insulating 
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FIG. 2: Conductivity a = Gx L/W and Fano factor F at the 
Dirac point (T^ate = 0), as a function of the aspect ratio of the 
graphene strip. The curves are calculated from Eq. Q in the 
limit N oo, for two different boundary conditions: smooth 
edge [solid curves, using Eq. ||1J] and "metallic armchair" edge 
[dashed curves, using Eq. 0]. The limit W/L oo (dotted 
lines) is given by Eq. @, regardless of the boundary condi- 
tion. The data points for the armchair edge are the result of a 
numerical solution of the tight-binding model on a hexagonal 
lattice. 



to metallic in the limit W/L but has no effect 
in the opposite limit W/L ^ oo, cf. Fig. |21 

To test the analytical results from the continuum de- 
scription of graphene, we have also carried out numerical 
simulations using the tight-binding model with nearest- 
neighbor hopping on a honeycomb lattice with metallic 
armchair edges. We took some 2000 lattice sites for the 
graphene strip, coupling it to semi-infinite leads at the 
two ends. The valley degeneracy of modes n = 1,2,... 
is now only approximate, but the relative magnitude of 
the mode splitting vanishes cx a/W as the width becomes 
large compared to the lattice spacing a |l7j . The numer- 
ical results, included in Fig.|21 are in excellent agreement 
with the analytical prediction. 

Fig. 121 shows that the minimum in the conductivity 
at the Dirac point is associated with a maximum in the 
Fano factor. The limiting behavior at the Dirac point for 
a short and wide strip is 
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cr^go/TT, F->l/3, for W/L ^oo. (9) 
Note that these limits are already reached for moderate 
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FIG. 3: Gate voltage dependence of the conductivity and 
Fano factor for a fixed aspect ratio. The conductivity min- 
imum at the Dirac point corresponds to maximal Fano fac- 
tor. The curves are calculated from Eq. for the case of 
a smooth edge (solid curves, taking q„ = (n -I- l/2)n/W) or 
metallic armchair edge (dashed curves, taking q„ = mr/W). 
The oscillations signal the appearance of propagating modes 
in the graphene strip with increasing gate voltage. 



aspect ratios W/L > 4. We derived this limiting behavior 
for two types of boundary conditions (smooth edge and 
metallic armchair edge), but we are confident that the 
result is universal, in the sense that it holds for the most 
general boundary condition at the edges of the graphene 
strip (as classified in Ref. Tg'I). 

The result ((SJ for the minimal conductivity agrees with 
other calculations |M II, which start from an un- 
bounded disordered system and then take the limit of 
infinite mean free path I. There is no geometry depen- 
dence if the limits are taken in that order. In the ballistic 
bounded system considered here, a geometry dependence 
persists in the thermodynamic limit |^. Existing exper- 
iments 1, are quasi-ballistic, with I ^ W < L, finding 
a ^ go- This would be consistent with Fig. [5^ if the ef- 
fective length in the experiment is set by I rather than L. 
Residual disorder may also explain why the oscillations 
in Fig.|21i, associated with the appearance of propagating 
modes, are not observed in the experimental gate voltage 
dependence (compare figure 1 of Ref. 0). 

The limit = 1/3 for the Fano factor is smaller than 
the value F — 1 expected for a Poisson process. The same 
1/3 value appears in a disordered metal 0,0, where it is 
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a consequence of classical diffusive dynamics. This corre- 
spondence is remarkable, since in our ideal graphene strip 
the classical dynamics is ballistic. The relativistic quan- 
tum dynamics of confined Dirac fcrmions is known to 
exhibit a jittering motion called "Zittcrbewegung" [Tsjl . 
originating from the interference of states with positive 
and negative energy 0, ^| . Our calculation implies that 
this relativistic quantum dynamics produces the same 
shot noise as classical diffusion. 

In conclusion, we predict that electrical conduction 
through an ideal graphene strip is associated with time- 
dependent current fluctuations — at zero temperature 
and without any impurities or lattice defects. The elec- 
trical noise is largest, relative to the mean current, when 
the Fermi energy is adjusted such that electrons and holes 
are degenerate. At this Dirac point the Fano factor (ra- 
tio of noise power and mean current) takes on the uni- 
versal value 1/3 for short and wide strips. Observation 
of this quantum-limited shot noise would be a unique 
demonstration of electrical noise produced by relativistic 
quantum dynamics. 
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lowship from the Polish Science Foundation. 



APPENDIX A: CALCULATION OF THE 
TRANSMISSION PROBABILITY 

1. Mode matching 

Plane wave solutions of the Dirac equation 
hv f(Jxdx + (Jydy 







take the form 



dx (^y 



g 1 ■9 + eV9 = e^ (Al) 



Xn,k{y) = a„ I I e"?"^+a' 



(A2) 



b„ \ J I e-*«"^ + bl 



The complex number Zn^k is defined by 

k + iqn 




Zn,k i 



(A3) 



(A4) 



where the + sign applies to the conduction band and the 
— sign to the valence band. For later use we note the 
identity 

Zn,kZn-k = -1- (A5) 

The energy of the state in the leads (x < and x > L, 
where V — Voo) is given by 



e^eVoo± fiv^k^ + ql (A6) 

In the strip (0 < a; < L, where V = Vgato) the wave 
vector k is replaced by fc, satisfying 



fc2 



(A7) 



(The ± sign again refers to conduction and valence 
bands.) In the leads we have only propagating modes 
(real k), while the modes in the strip may be propagat- 
ing (real k) or evanescent (imaginary k). The sign of k 
is not fixed by Eq. (|A7() and may be chosen arbitrarily. 

The transverse wave vector g„ as well as the coeffi- 
cients a„, a^, 6„, b'n of the n-th mode are determined (up 
to a normalization constant) by the boundary conditions 
at y = and y = W. We consider a class of boundary 
conditions for which the resulting parameters are inde- 
pendent of the longitudinal wave vectors fc, k. Because 
Qn is then the same in the strip and in the leads, the 
mode matching at the lead-strip interface does not mix 
the modes and can be done analytically. 

A scattering state incident from the left at energy e 
in the n-th mode is constructed from the states (IA2I). 



(A8) 



<i>L if a; < 0, 

ifO<a;<L, 
$K if X > L, 

= Xn.fce''^"' + r„XrK-/ce"''='', (A9) 

<i>fl = te,fce''("-^\ (AlO) 

^ = «nX„,^e»'"+/3"X„,_^e-^^ (All) 

The sign of the (real) wave vector k is chosen positive for 
the conduction band and negative for the valence band. 

The reflection and transmission amplitudes r„ , i„ and 
the two coefficients q;„,/3„ are determined by requiring 
continuity of at x = and x — L. There are four 
independent equations. 



-I3n 



1 

'n,-k 



(A12) 



-ikL 



(A13) 



Notice that the coefficients a„ , aj^ , 6„ , b'^^ have dropped 
out. We need the solution for the transmission ampli- 
tude. 



(i + 4fc)(i + <fc) 

e'^^(z„,fc - z^ ^)2 + e-^'^^il + Zn,kz,^ },) 



(A14) 
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where we have used the identity (|A5|) . The transmission 
probabiUty r„ of the n-the mode at the Fermi level is 
obtained from r„ — |i„p at e = 0. 

We are interested in the limit |Voo| — > oo of an infinite 
number of propagating modes in the leads. Then Zn^k — > 
1, hence 



with (/) = 0, one finds 



mr an _ Q,i _ -. 



(A21) 



2 + 2z" 



k cos(fcL) + «(eV^ato/?iw) sin(/cL) 



(A15) 



with fc^ = (eVgato/^w)^ — (Z^i- This result corresponds to 
Eq. jSJl in the main text. 



2. Boundary conditions 

It remains to calculate the transverse wave vector g„, 
for specific boundary conditions at y = and y = W. 
Substitution into Eq. ljA15(l then completes the calcula- 
tion of the transmission probabilities. 

For the infinite mass boundary condition, 



^|,=o= (''n^ _l )*|,=o, (A16) 



-a, 

cr 



^y^w = ( a ) my^w, (A17) 



The modes with n = 1, 2, 3 . . . have a twofold degeneracy 
(because the choices a„ = and a'^ — produce two 
independent states). The lowest mode n — 0, however, is 
nondegenerate (because l^o^fcl = 1, so the states produced 
by a„ = and = are the same). 

A semiconducting armchair edge (corresponding to a 
strip width which is not an integer multiple of three unit 
cells) has = ±2tt/3 in Eq. (|X20|l . Then aU modes 
are nondegenerate, consisting of one sequence with q„ = 
(n — (f>/2n)Tr/W, a„ = 5'^ = 0, = bn and one sequence 
with q„ = {n + (f)/2T:)7:/W , a'„ = bn = 0, a„ = b'„. 

The results for the conductance and Fano factor at the 
Dirac point, including all degeneracies, may be summa- 
rized by 



2p2 



E r„(i-T„) 



n— — oo 



E Tn 



-,(A22) 



n— — OO 



1 



cosh^[(n-|-a)7rL/Vl^]' 



(A23) 



one finds 

Each mode n = 0,1,2,... has a twofold valley degener- 
acy. 

For the metallic armchair edge. 



(A19) 



e*"^ 




^\y^w, (A20) 



with a = 1/2 for the infinite mass boundary condition, 
a = for the metallic armchair edge, and a = 1/3 for 
the semiconducting armchair edge. 

The general formula (|A23p holds under the assumption 
that longitudinal and transverse momenta are not cou- 
pled by the boundary condition. The zigzag edge is an 
example of a boundary which does couple k and g„. In 
that case the jump in fc at x = and x — L introduces a 
mode mixing which complicates the analytical solution. 
We do not expect any differences to appear in the limit 
® for W/L oo. 
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